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1. INTRODUCTION 

The performance of high data rate transmissions over wireless fading channels severely degraded due 
to the multi path effects which causes inter symbol interference (ISI). In order to combat the fading effects, 
OFDM has been widely adopted to wireless transmission [1-8]. The multipath propagation causes a time 
varying channel state information (CSI), which needed to be predicated or estimated using channel estimation 
techniques in order to recover the transmitted signal. Pilot aided channel estimation is the widely used 
technique begins from the traditional techniques such as, least square (LS) and linear minimum mean square 
(LMMS), ending with many recent ones used to improve the estimation performance [2, 9]. 

Compressive sensing (CS) is one of the recent techniques adopted for channel estimation in OFDM 
systems by exploiting channel sparsely representation with dictionary basis [10]. Several approaches have been 
employed to construct matrices in order to represent the channel in a sparse manner such as, discrete Fourier 
transform (DFT) and random dictionaries. These approaches do not consider the time variation property of 
the channel since the time varying channel parameters aren't taken into account [11-14]. 

The main contribution of this paper is, the design of a novel delay dictionary-based CS technique to 
overcome the problem of the dictionary proposed by [15] in estimating LTV channel in the presence of Doppler 
frequency shifts. In [15], a sample spaced delay dictionary was proposed to recover the CSI using CS in 
multiple input multiple output (MIMO)-OFDM system. The concept of the research based on estimate 
the channel coefficients for a time varying channel considering the useful OFDM symbol duration regarding 
the guard band, and tacking the delay profile into account. Considering channel delay parameters, the dictionary 
proposed by [15] improves its ability to recover the CSI even when number of pilots reduced, while it fails 
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in estimating LTV channel coefficients in the presence of Doppler effect since it exploits the time variability 
characteristic of the channel. Which leads to a conclusion that, the dictionary proposed by [15] can't be applied 
for channel estimation in LTV channels since it doesn't sense the Doppler effect of the channel. 

The rest of this paper is organized as follows, in section 2, a brief introduction to CS theory introduced 
with the required analysis of LTV channel and the proposed system model for sparse channel estimation. 
Simulation tests, required system parameters, and test results are introduced in section 3. Finally, the main 
concluded remarks and future work are listed in section 4. 


2. LTV CHANNEL ANALYSIS AND ESTIMATION BASED CS THEORY 
2.1. Compressive sensing 

Since the idea behind signal sparsity appears, many publications of sparse signal representations and 
compressive sensing introduced especially in signal processing community [16]. With compressive sensing, 
areal finite signal x € R™, can be expressed in an orthonormal basis; 


x = 21 yib: (1) 


where Y = [YY ... Wy] represents the orthonormal basis, and 0 = [6,02 ... Om| is the sparse vector where 
the number of non-zero elements (K<<M) much smaller than the number of zero elements and named as 
a K-sparse vector. Using matrix notations, x = W@, where Y of size M x M [17]. Consider a classical linear 
measurement model where y = x = WO. Where O represent the k-sparse vector of size M x 1 to be 
estimated using the effective measurement matrix dyp, where the measurement matrix @ is of size N x M, and 
y is the measurement vector of size N x 1. Hence, each observation of y vector represents the projection of 
vector x on a row of the sensing matrix ¢ as described in Figure 1 [18]. 

From the mathematical expression of CS in the Figure 1, it is clear that a non-linear system of 
equations must be solved to recover the sparse vector 0, where the number of observations N is much less than 
number of unknowns M. Since ¢@ matrix projecting the vector x, low value of incoherence is required to insure 
mutually independent matrices and hence better CS performance. The maximum value amongst inner product 
of the Orthonormal basis and the orthonormal measurement matrix defined as incoherence. Therefore, to 
recover the sparse vector correctly from y = @x, the sensing matrix @ should be designed carefully [19, 20]. 





: m Mx 1 
i = sparse 
signal 
K 
NxI —_— -— » nonzero 
measurements | | entries 


y = Yo 


Figure 1. CS mathematical representation 


2.2. System model and LTV channel 

In this paper, the OFDM system of Figure 2 is considered. At the transmitter Ty side of this system, 
a stream of symbols x[k] (data d[k] and pilots p[k]) are mapped using binary phase shift keying (BPSK), where 
x[k] split into data blocks after serial to parallel conversion. Each of these blocks represent OFDM block 
contain data and pilot symbols. The length of each OFDM block is N subcarriers. A cyclic prefix (CP) of length 
(Lep) is prepended to each OFDM block to prevent adjacent interference and considered as a guard band (Gj). 
After CP insertion, the OFDM block transmitted over an LTV channel which is a multipath propagation 
channel. In the proposed work, the LTV channel has been assumed to have a finite impulse response with L 
paths. The transmitter and receiver are assumed synchronized in both time and carrier frequency. 

The multipath fading channel response is expressed as follows [21, 22]; 


h(t) = Xio a * O(t — Ti) (2) 
where, the ith path of wireless environment is characterized by a propagation delay (t;) and attenuation 
(ai). The received baseband signal r(t) is modeled by two components, amplitude and phase, which can be 


expressed as; 
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r(t) = Xit a; * E(t — Ti) sere (3) 


where, e~727fcTi is the complex phase factor, and for a narrow band transmission; 
(t — Ti) = (t) [23]. Hence; 


h= gio are AT (4) 
The time varying property of the channel implies that channel coefficients changed over time. This changing 


is related to the change in the frequency of the received signal, which related to the relative movement between 
the transmitter and the receiver, Hence, the corresponding channel delay (t;) is changing [21, 22]; 


T(t) = q -1 (5) 


where T;(t) is a function of distance, and hence; 


vcos 6t 





A(t) = a di o eel (6) 
— ŞL-1 p. p-j2īfcTi jenje SR 
h(t) 5 aj e e c (7) 
where, f, —— ” is the Doppler frequenc = cos 0, where = f, is the maximum Doppler shift 
c pp q Y Ja dmax dmax Co Pp 


(famax). Assuming that the movement of the mobile system is uniformly distributed from 0 < @ < m rad, and 0 
is normalized. Thus, the channel impulse response 1s; 


A(t) = DEA a; eriet elmrat (8) 


Since this component (e/2"/4*) is a function of time, as a result, the channel coefficients A(t) are time varying. 
Such a time varying channel is known as a time selective channel. How fast or slow the channel changes 

depends on the channel coherence time (T) where the channel is approximately constant during 
T, [21]. 

At the other hand, in order to estimate the time variant channel coefficients using CS technique, 
the sensing matrix should be designed with atoms related to the two effecting parameters (t; and fa) of 
the LTV channel of (8). This will lead to compute the rate of change of the wireless channel by analyzing 
the correlation between channel coefficients. Assume that a;(t) is the channel coefficient at the ith path 
at time t; 


a; (t) = A; e IJ2ufT; el2tfat (9) 


Thus, to compute the correlation between a;(t) and a;(t + At), the expectation of a;(t) and a,(t + At), 
E{a;(t) a;(t + At)} should be computed [24]; 


a;(t) = a; e J2mfti ei2mfat (10) 

a;(t + At) = a; e~J27FT% es2mfalttAt) (11) 
thus; 

W (At) = E{|a;|? * e/27fa44 (12) 


where ¥ (At) refer to the correlation function between a;(t) and a;(t + At). Let |a;|* normalized to be 1, thus, 
W(At) = Efel*™faat} (13) 
which summarized as; 


P (At) = Jo(2mfaAt) (14) 
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where f4 represents the maximum Doppler frequency and Jọ is the Bessel function of 0“ order. Finally, 
the autocorrelation function ¥ (At) of LTV channel can be expressed in terms of coherence time [25]; 


V(At) =G. (15) 


where; At = xT,,x = 1,2,3 ....... , and, T, =— 


s[k] 
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Figure 2. OFDM System Model 


CP removing 


2.3. Sparse channel estimation 

Since CS has gained a much popularity in communications, recently, it is one of the small numbers of 
strong paths used for channel estimation. It assumes that sparse signals can be approximated with a small 
number of measurements compared to the large number required with Shannon-Nyquist rate [1]. Hence, to 
estimate the channel vector h € C™*1 from y measurements, a CS problem of section (2.1) should be solved, 
where y is expressed as follows; 


y=Ah+n (16) 


N 
n: AWGN noise with zero mean and variance a, = = 


A: The sensing matrix 

The sparse representation of data in terms of atoms is the main objective of the dictionary design, 
which later used to reconstruct the sparse signal, where h assumed to be K-sparse CSI and its energy uniformly 
distributed among a small number of taps without any prior knowledge of their location, which must be 
estimated with effective sensing. It is clear from (15) that channel coefficients are changed with respect to the 
channel coherence time. Therefore, in this paper, the dictionary matrix is designed in a manner in which the 
two delay parameters of the autocorrelation function are taken into account. The equispaced pilot subcarriers 
p[k] are embedded within the data subcarriers d[k] of the OFDM system of Figure 2, where the number of 
training pilots is N,, and At is assumed to be a taped delay profile along the OFDM symbol. Where; 


At =|0,ix£,a] i=1,2,3,..,N (17) 
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GixTs 
ToXN 
which assumed to be the CP appended to each OFDM symbol in order to mitigate the ICI, and T; is the OFDM 
sample time. An N x N dictionary matrix is constructed with atoms related to each subcarrier position £; along 
the OFDM block length Tsymbo and multiplied by the taped delay atoms qt; of At. 
Tsympol = (N + Gi) Ts, is the OFDM symbol time including G;. Therefore, the dictionary Dyxy is represented 
as follows; 





a represent the minimum channel tap spacing which equals to (G; xX T; — ), and G; is the guard interval 








_j2n £474 _j2n €4TN 
e Tsymbol eis e Tsymbol 
D=| : ; (18) 
_j2metnty _J2T NTN 
T T 
e Symbol soe e ‘Symbol NXN 


where the rows of dictionary matrix D are refer to the subcarriers positions along the OFDM symbol, while 
columns are refer to the delay vector of each subcarrier. Regarding sensing matrix. A construction, an Np rows 
are selected from D related to pilot locations, and multiplied by N, X N matrix of pilot data using dot 
product multiplication; 











_ j2m e477 _j2mt1tN 
e Symbol ... e TSymbol eI Py ees e It P1(N) 
A= : k | i: (19) 
_J2T tN D1 _J2T NTN e IT PNp e I” PNp(N) NXN 
e TSymbol s. -ø Tsymbol p 


NpXN 


P; is the pilot symbol used for estimation considering equally likely symbols of [0, 1], and 
i = 1,2,3, ...., Np. Different sparse signal recovery algorithms can be applied to solve the CS problem in 
a number of iterations to minimize CS error with respect to D. Once the channel dominant taps estimated 
(i.e. recovering of the sparse vector 0), the whole CIR is built at all locations simply from h = D x 0, and 
the equalization process implemented. 


3. SIMULATION TEST AND RESULTS 

In this paper, the OFDM system of Figure 2 compares the test results of both least square (LS) and basis 
pursuit (BP) based channel estimation techniques. Different OFDM system parameters are listed in 
Table 1. In addition to AWGN noise, a 6 tap LTV channel is considered as a Rayleigh fading channel with paths 
delays and power vectors standardized by ITU channel model of Table 2. Basis pursuit (BP) algorithm is used to 
solve the convex optimization problem with MatLab for sparse signal recovery, where it uses the l,_ norm to 
regularize the problem [9]. The performance test of OFDM system was shown in the form of bit error rate (BER) 


versus signal to noise ratio (SNR), where SNR is determined by the corresponding (Fo /N ) in dB. 


Table 1. OFDM system parameters 


Parameter Value 
Number of transmitted bits 64000, 128000, 256000 
Modulation BPSK 
Sampling Time (T;) | usec 
OFDM Subcarriers 64, 128, 256 
Number of pilots 16, 13 
Cyclic prefix length (Lep) 16 


Maximum Doppler shift (fz)Hz 9, 10, 20, 40 


In Figure 3, the BER performance of OFDM system over a multipath indoor channel environment 
using LS and BP is presented. In addition to compare the estimation performance of BP over LS algorithm, 
the purpose of this test is to show the ability of the proposed dictionary to recover the CSI with different 
delay parameters for both indoor and outdoor environments. In the present test, BP performance outperforms 
LS technique by about 4.5 dB at a BER of 107° with 16 pilots out of 64 subcarriers and zero Doppler 
frequency. In Figure 4, LS and BP algorithms are tested over a multipath outdoor channel environment of 
Table 2. In this test, the OFDM block containing 64 subcarriers and the number of pilot subcarriers used is 
16. Different Doppler shifts are considered in order to test the recovering ability of the proposed dictionary 
in the presence of Doppler effects. 
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Figure 3. BER performance with Figure 4. BER performance with 
fa = 0 Hz, Np = 16,and N = 64 N, = 16,and N = 64 


Table 2. ITU Channel Models [26] 


Indoor Outdoor 
Delay (ns) Power (dB) Delay (ns) Power (dB) 

0 0 0 0 

50 -3 310 -1.5 
110 -10 710 -9.0 
170 -18 1090 -10.0 
290 -26 1730 -15.0 
310 -32 2510 -20.0 


As could be noticed, CS based channel estimation algorithm improves the estimation performance 
as compared to LS algorithm even with Doppler effect. By comparing BP performance for both fg = 0 Hz, 
and f4 = 10 Hz, it is clear that as fg increased to 10 Hz, the performance test degraded by about 15 dB at 
a BER of 107%. At the other hand, BP performance degraded when f4 increased more than 10 Hz and become 
worse than LS unless the number of pilots used for estimation increased, which in turn improves 
LS performance. The same test was repeated with a lower number of pilots, where 13 pilots was inserted within 
the OFDM block at equally spaced locations instead of 16 pilots as shown in Figure 5. The test shows that 
the BP performance degraded but still much better than LS. This observation leads to the possibility of using 
reduced number of pilots for channel estimation without sacrificing the accuracy of channel estimation, when 
the rate of change of channel coefficients increasing according to Doppler shift effects. 





0 5 10 15 20 25 30 35 40 


Figure 5. BER performance with N, = 13,and N = 64 
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Another test considering different subcarrier numbers and Doppler shifts is shown in Figure 6. 
The test results proved that BP exceeded LS performance. But this superiority is still limited by the amount of 
Doppler shift and number of pilots, where it is degraded when f4 increased above than 20 H, and 40 H, for 
128 and 256 OFDM subcarriers respectively with 16 pilots. This degradation is shown in Figure 7, where LS 
and BP algorithms are tested with, fg = 30H, for N = 128, and f4 = 50H, for N = 256. Finally, 
it could be concluded that, as the amount of Doppler shift increase, the estimation performance degraded due 
to the Doppler effect on the channel. This degradation manifests itself when the number of subcarriers 
increased, where the subcarrier bandwidth will be decreased, so, it is more sensitive to the Doppler and requires 
an additional process to eliminate carrier frequency offset (CFO). 


—+— BP-N=128,f =20 | 
—a— LS-N=128,f,=20 
—e— BP-N=256,f =40 | 

—#— LS-N=256,f,=40 
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Figure 6. BER performance with Figure 7. BER performance with N, = 16 
N, = 16, N = 128 and 256 


4. CONCLUSIONS AND FUTURE WORKS 

In this paper, the proposed dictionary design was tested to achieving the desired results of BP based 
CS algorithm in estimating of CIR of a LTV channel. At the other hand, this performance is limited to the low 
to moderate Doppler frequency shifts. The future work may be carried to extend the current work to be used 
for estimation of LTV channel with high mobility or high Doppler frequency shifts. 
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